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The source of spin-wave resonance (SWR) in thin films of the ferromagnetic semiconductor (Ga,Mn)As is 
still under debate: does SWR stem from the surface anisotropy (in which case the surface inhomogeneity (SI) 
model would apply), or does it originate in the bulk inhomogeneity of the magnetic structure of the sample 
(and thus requires the use of the volume inhomogeneity (VI) model)? This paper outlines the ground on 
which the controversy arose and shows why in different conditions a resonance sample may meet the 
assumptions of either the SI or the VI model. 

Dilute ferromagnetic semiconductors are a class of very promising materials of the future 16 . Gallium 
manganese arsenide (Ga,Mn)As, created on the basis of the semiconductor gallium arsenide by the 
addition of a small percentage of manganese as a magnetic dopant, is one of the most intensively studied 
compounds in this class 715 . The free motion of positive charge carriers (holes) throughout the crystal results in 
the ferromagnetic order of the manganese ions. The basic magnetic characteristics of the material depend on the 
amount of the dopant ions and the spatial distribution of the concentration of the charge carriers (holes) 
transmitting magnetic information between the Mn ions. A particularly interesting situation occurs in thin films, 
in which magnetic characteristics (magnetic anisotropy, magnetization, exchange length and stiffness constant, 
damping constant, etc.) are in general nonuniform along the growth direction perpendicular to the film surface. 
The character of this nonuniformity reflects the distribution profile of the charge carrier concentration in the film. 

The spatial magnetic profiles in thin films can be determined by means of ferromagnetic resonance, which 
reveals its fine structure in a multi-peak resonance spectrum in thin-film systems; this type of ferromagnetic 
resonance is referred to as spin-wave resonance (SWR), as each peak in the resonance spectrum corresponds to the 
excitation of a specific spin wave. On the other hand, the spectrum of allowed spin-wave excitations is determined 
by the shape of the magnon potential of the system. Since the position of each SWR peak corresponds to a spin- 
wave energy level resulting from the prevailing magnon potential, an experimental SWR spectrum can be turned 
into the corresponding profile of the magnon potential by an appropriate calculation procedure. Thus, providing 
information on the spatial distribution of the basic magnetic characteristics, including the charge carrier con- 
centration in the film, resonance measurements are of vital importance for the elucidation of the origins of 
ferromagnetism in the material under investigation. 

Spin-wave resonance in thin films has been studied particularly intensively in gallium manganese arsenide in 
the past decade 1216 27 . Especially rich resonance spectra were obtained in studies with a variable configuration of 
the static field with respect to the film surface. The field was rotated both perpendicularly to the film surface 
(which corresponds to variable polar angle 0 H between the direction of the external field and the surface normal) 
and in the plane of the film (variable azimuth angle tf> H between the external field and a reference direction in the 
film plane). The results of these measurements clearly indicate that the evolution of the SWR spectrum with the 
field configuration is correlated with that of the spatial distribution of the spontaneous magnetization and the 
anisotropy; thus, configuration and space dependence of the magnon potential should be assumed as well. 

In the present paper we shall only analyze SWR measurement data concerning the out-of-plane rotation of the 
magnetic field, mainly because of the controversy that arose in the interpretation of these results over an issue 
which therefore requires elucidation (in a separate paper we intend to analyze measurement data obtained in 
SWR studies with in-plane rotation of the magnetic field as well). If researchers tend to agree on the interpretation 
of SWR spectra in two extreme configurations - the perpendicular and parallel configurations, corresponding to 
0 H = 0 and 0 H = 90°, respectively - the interpretation of results obtained in intermediate configurations is under 
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debate. Almost as a rule, a particular configuration of the external 
field tends to occur in this range at a critical angle 0 H , for which the 
multi-peak SWR spectrum collapses to a single-peak FMR spectrum. 
There are two schools of thought regarding the interpretation of the 
occurrence of this critical angle. These two prevalent opinions agree 
on the physical state of the thin film in the critical configuration, but 
differ in the interpretation of the configuration- related processes that 
accompany the rotation. Both schools agree that in the critical con- 
figuration the thin film (its magnon potential, to be precise) is mag- 
netically homogeneous, and the boundary conditions (specifically, 
the surface spin pinning) correspond to the natural conditions, only 
resulting from the reduced neighborhood of the surface spins (a 
precise definition of the natural pinning conditions is provided in 
the next Section). The difference of opinion concerns the configura- 
tion evolution leading to the above-described "naturally homogen- 
eous" magnetic state. One school 21 uses the surface inhomogeneity 
(SI) model and assumes that rotation of the magnetic field does not 
modify the profile of the bulk magnon potential, which remains 
homogeneous across the film; only the surface pinning conditions 
change, diverging from the natural conditions as the angle grows 
above or decreases below the critical configuration (with the surface 
pinning decreasing or increasing). In contrast, the other school 27 , 
using the volume inhomogeneity (VI) model, claims that it is the 
bulk profile of the magnon potential that changes with the field 
configuration: remaining linear, but inclined at different angles with 
respect to the surface of the film, the magnon potential increases or 
decreases inside the film as the configuration diverges from the crit- 
ical angle, while the natural conditions prevail invariably on the 
surface. In this paper we opt for the interpretation based on the SI 
model and propose a theoretical model of the configuration evolu- 
tion of the surface spin pinning in agreement with the experimental 
data. Our interpretation leads to some physical conclusions, which 
provide new insights into the surface properties of ferromagnetic 
semiconductor (Ga,Mn)As thin films. 

The Goal of the Study and the Concept of SWR 
Surface Pinning Parameter 

Our discussion of the state of the art of the research in the critical 
angle effect in SWR in ferromagnetic semiconductor (Ga,Mn) As thin 
films will rely on the representative study performed by Liu et al, 
reported in Ref. [21], presenting SWR spectra measured for inter- 
mediate angles 0 H between the external field and the surface normal. 
Characteristically, in the out-of-plane configuration, with the field 
vector rotated in a plane perpendicular to the surface, the SWR 
spectrum, consisting of multiple peaks in the perpendicular (0 H = 
0) and parallel (0 H = 90°) configurations, is found to collapse to a 
single-peak FMR spectrum in an intermediate configuration corres- 
ponding to a critical angle 0 C H (19° in the studied sample). 

The critical angle effect in SWR has been known for years, but that 
observed in (Ga,Mn)As samples is very unusual. The peculiarity is 
that the critical angle 0 C H coincides with the border between two 
configuration domains in which the SWR spectrum fulfills the 
assumptions of different models: the surface inhomogeneity model 28 
for Oh > 0 C H (in which range the spacing between the resonance peaks 
is proportional to n 2 , where n is the spin-wave mode number), and 
the volume inhomogeneity model 29 , which applies for 6h < 0 C H 
(where the spacing between the resonance modes is proportional 
to n). A question arises: what mechanism underlies the occurrence 
of the inhomogeneity, if surface inhomogeneity prevails for 0h>0 c h , 
and volume inhomogeneity for 0h<0 h 7 . And what particular surface 
mechanism leads to the occurrence of the critical angle 0 H at which 
these two types of inhomogeneity fail to be "seen" in the resonance? 
In the present paper we are going to find an answer to this question 
by analyzing in detail the surface anisotropy fields found in 
(Ga,Mn)As thin films in studies of SWR in this material. 



The experimental SWR spectra analyzed in this paper were mea- 
sured in the "out-of-plane geometry", as referred to by the Authors of 
Ref. [21]. In this out-of-plane geometry, the (Ga,Mn)As layer was 
cemented to a parallelepiped of GaAs (100) substrate material, the 
[110] edge of the specimen oriented vertically. The external magnetic 
field Hwas confined to the horizontal plane (i.e. perpendicular to the 
film surface) allowing SWR measurements with Hin any intermedi- 
ate orientation between the normal to the film surface, H\ \ [001], and 
the in-plane orientation, ff||[ll0]. In this particular geometry the 
magnetization M of the sample lies in the same horizontal plane as 
the field H. Thus, the spatial orientation of the vectors H and M is 
defined by two polar angles, 0 H and 0 M , between the respective 
vectors and the normal to the surface of the film. For (Ga,Mn)As 
samples in this particular geometry of the external field the free 
energy density F ± of the system has the form 24 : 

F± = —M x [ — 2H (cos Om cos Oh + sin Om sin Oh ) 

+ (4rcM-H 2 i)cos 2 0 M — ^#4± cos 4 0 M (l) 
— -.tfyl sin 4 C M — H 2 \\ sin 2 0 M , 

where H 2± and H 4± are the uniaxial and cubic anisotropy fields, 
respectively, perpendicular to the plane of the sample; ff 2 || ana #411 
are the in-plane uniaxial and cubic anisotropy fields, respectively. 

The above equation includes as many as four anisotropy fields. 
Although epitaxial (Ga,Mn)As films crystallize in the cubic zinc- 
blende structure, which only justifies the presence of the two cubic 
fields, H 4 i and H 4± , in Eq. (1), the lattice mismatch with the sub- 
strate distorts this initial structure and reduces its symmetry from 
cubic to tetragonal, which results in the appearance of two more 
fields, H 2 \\ and H 2± , related to the uniaxial anisotropy. However, it 
should be also considered that the crystal symmetry is additionally 
broken at the surface; as a result, both the uniaxial and cubic terms 
differ by some amount in the bulk and surface regions (as confirmed 
by SWR experiments, see Ref. [30]). This implies resultant surface 
anisotropy fields essential for the surface spin pinning mechanism 
responsible for the spin dynamics features observed in (Ga,Mn)As 
films. Another important source of surface anisotropy in (Ga,Mn)As 
films is believed to be the fact that in this material the charge carrier 
(hole) concentration on the surface is different than in the bulk 10,3132 , 
which also affects the values of all the above-mentioned surface 
anisotropy fields. 

It should be noted that the SWR studies of (Ga,Mn)As conducted 
so far tended to focus on volume characteristics only, such as the 
uniaxial anisotropy or the exchange constant of the studied material. 
The aim of this paper is to use SWR for getting a better insight into 
the ferromagnetism of dilute semiconductors in terms of their sur- 
face characteristics, the current knowledge of which is scarce. For this 
reason, in the analysis presented in this paper, we refer to our earlier 
quantum theory of SWR 33 ~ 42 , in which we have introduced the con- 
cept of surface spin pinning parameter, a quantity that measures the 
degree of pinning of the surface spins and reveals explicitly different 
surface magnetic anisotropies present in thin films. 

The concept of surface pinning is related to the description of the 
energy status of surface spins, specifically to the degree of freedom of 
their precession. In a very simplified image introduced in Refs. [34, 
39], besides the effective magnetic field present throughout the sam- 
ple, an additional magnetic field K sur f, referred to as the effective 
surface anisotropy field, acts on the surface spins. As we have shown, 
the boundary conditions to be fulfilled by the precession of the sur- 
face spins can be expressed by the surface pinning parameter, 
defined: 
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A=\- 



rftlt, 



(2) 



where a is the lattice constant, D ex is the exchange stiffness constant, 
and m denotes a unit vector oriented along the magnetization M of 
the sample. Note that a complete lack of anisotropy field on the 
surface corresponds to the surface parameter value one; the freedom 
of the surface spins in this situation will be referred to as the natural 
freedom. In the case of nonzero anisotropy field three situations, 
substantially different from the physical point of view, may occur 
depending on the angle between the magnetization Mand the surface 
anisotropy field K sur f. If the surface spins are aligned perpendicularly 
to K SU rf, their freedom remains natural (A = 1); otherwise, the sur- 
face spins are pinned (and A < 1) or unpinned (and A > 1) for the 
above-mentioned angle acute or obtuse, respectively. All three pin- 
ning regimes are schematically depicted in Figure 1. In the rigorous 
theory of SWR the surface pinning parameter can be represented 
(Cracknell and Puszkarski 37,38 ) as a series expansion in spherical 
harmonics Yj m (0, <j>): 



A(9,<j>) = l- "-K surf (0,<l>)-m=J2 A 'mY, m (0,< 

LJ. - 



/ = 0 m=-l 



(3) 



1 = 0 



aiP°(cosO) P™ (cosO)x(a im coswj</> + f> im sin m<j>) 



1= -/ 



where 0 and (j> are the out-of-plane polar angle and the in-plane 
azimuth angle, respectively, of the magnetization M. The coefficients 
ai, 0!j m and fti m (which can be found experimentally) determine the 
respective energy contributions brought to the effective surface pin- 
ning by different surface interactions. As established in Ref. [38], in 
the case of surface cut (100) - which is that of the thin-film samples 
considered in Ref. [2 1 ] - all the terms with odd values of Z vanish, and 
the only values allowed to m are 0, 4, 8, .... In our research we have 
also observed 39 that in the case of thin films the series (3) can be cut to 
only include terms up to I = 4, since further contributions tend to be 
minor. Thus, we propose the following angular dependence of the 
surface parameter as appropriate for the interpretation of the SWR 
spectra obtained in Ref. [21]: 

A(8,<j>) = 1 -o 0 - a 2 {0,<t>) (3 cos 2 0- l) -a 4 (0^)cos4f (4) 

The above formula provides the basis for the elucidation of the most 
important surface mechanisms behind the SWR surface dynamics in 
(Ga,Mn)As thin films, which is the main goal of the present paper. 

Determination of the Equilibrium Direction of 
Magnetization in (Ga,Mn)As Thin Films 

In the investigated (Ga,Mn)As sample four bulk anisotropy para- 
meters appearing in the Eq. (1) have the values 21 : 



4nM eff = 4nM-H 21 =4588 Oe, 



= 0, Hm = 197 Oe, H 2 



= 77 Oe. 



(5) 
(6) 



Let us determine now the equilibrium direction of the magnetization 
of the sample, i.e. the equilibrium angle 0 M . We will use the condition 
of equilibrium of the system, which requires the first derivative of its 
free energy F ± to vanish: 

dFj_ 



= 0; 



(?) 



this condition allows to determine the sought relation 0 M = 0 M (Oh)- 
Since the condition (1) must be fulfilled when resonance occurs, 
the magnetic field H in ( 1 ) is the resonance field, H = H res ; we read its 
value from Figure 5 in Ref. [21], identifying it with the resonance field 




b) A< 1 







w \ M 
■^surf \ / 
















FILM 






Figure 1 | Schematic representation of three surface spin pinning 
regimes which prevail in a thin film depending on the configuration of its 
magnetization M with respect to the effective surface anisotropy field 
Ksurf ( see (2)). When aligned as in (a), the surface spins do not feel the 
anisotropy field and A = 1, which corresponds to their natural freedom. In 
the configurations (b) and (c) the surface spins are pinned (A < 1) and 
unpinned (A > 1), respectively, due to the anisotropy field. 



of the fundamental mode (n = 1). The 0 M = 0 M (0 H ) relation deter- 
mined numerically on the basis of the above considerations is shown 
in Figure 2; we refer to this relation many times in this paper when 
analyzing the experimental SWR spectra reported by Liu et al. 21 . 

Surface Versus Bulk Uniaxial Anisotropy 

In this section we shall consider the case in which only the perpen- 
dicular uniaxial anisotropy H 2 ± enters the formula (1) for the free 
energy. In that case the free energy reads: 
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Fx = -M x [ — 2H(cos6 M cos 0 H + sin 0 M sin9 H ) +4nM e g cos 2 0m] , 
and the use of the well-known Smit-Beljers resonance formula: 



1 



M 2 sin 2 0 M 



8 2 F, 8 2 F, 



d 2 F, 



8(j> M dO M 



(9) 



leads to the following configuration resonance condition, only apply- 
ing to the uniform mode k ± = 0 in the case considered: 

2 

= [Hcos(9 M — 6 h )— 4izM e ffCos20 M ] ^ 
X [H cos(O m - 0 H ) - $nM e fj cos 2 0 M ] ■ 



It will be very informative to derive the same condition in the micro- 
scopic approach, in which the energy of the system is expressed by 
the Hamiltonian: 

ij^tf ij ij 

its successive terms account for the isotropic exchange interaction, 
the Zeeman energy of the spins, and the perpendicular uniaxial 
anisotropy energy. The subscript Ij defines the position of the given 
spin, with I labeling the layer and the two-dimensional vector j defin- 
ing the position of the spin Sy in the l-th layer. The energy of a 
standing spin wave with a wave number k ± in this model is given 
by the expression 44 : 



2DS 

Hcos(0 M — Oh) H cos20 M ■ 

g^B 



2Sz±Ja 2 

g^B 



Hcos(0m — 6h) ■ 



IDS "'-■'" 2 
cos 0m + 

g^B g^B 



(12) 



2Sz±Jd 2 



This condition is the counterpart of the condition (10) obtained in 
the macroscopic approach (for k ± 0). From the comparison of 
these two formulas it follows that: 



4nM, 



2DS 



4i = 



and the coefficient at k\ can be identified as: 



D ex = 



2Sz±Ja 2 

g^B 



(13) 



(14) 



CD 

c 

£= 
o 



01 



=5 

cr 




10 20 30 40 50 60 70 80 90 
Magnetic field angle, 8 H 

Figure 2 | Equilibrium magnetization angle 0 M vs. the external field angle 
ff H as determined from the condition (7) for the (Ga,Mn)As thin film 
studied by Liu et a/. 21 ; 9 C M and 8 C H are the respective critical SWR angles. 



To obtain the formula for the surface parameter expressed by the 
surface perpendicular uniaxial anisotropy, we must yet rewrite the 
third term in the Hamiltonian (11) in the generalized form: 



where the uniaxial anisotropy constant D; is assumed to be: 
D s for surface spins, 



Db for bulk spins. 



(15) 



(16) 



On the basis of our earlier papers 42,45 it can be demonstrated that in 
the approximation assuming circular spin precession the following 
expression for the surface parameter results from this model: 



D b -D s , ,„ 
A = l-^— - - l-3cos 2 0 A 
2zJ v 



(17) 



where D b and D s , as indicated above, denote the bulk and surface 
values, respectively, of the microscopic uniaxial anisotropy constant. 
Now, using the identity relations (13, 14) we obtain the sought final 
formula in which the surface parameter is expressed by macroscopic 
quantities: 



A(0 M ) = 1- 



1 



4ti M 



-surface 



(l-3cos 2 tf M ) 



(18) 



Out-of-Plane Angle Dependence of the Surface 
Parameter in (Ga,Mn)As Thin Films 

In the present paper we shall focus on the configuration dependence 
of the SWR spectrum of (Ga,Mn) As thin films with the external field 
ffonly rotating in a plane perpendicular to the surface of the sample 
from the direction along the surface normal (0 H = 0) to the in-plane 
direction {0 H = 90°). According to the formula (4), in this case the 
surface parameter of a (Ga,Mn) As thin film can be represented as the 



A(6 M ) = 1 - flo - a 2 (0 M ) (3 cos 2 0 M — l), 



(19) 



where 0 M is the angle between the surface normal and the magnet- 
ization M of the film (let us remark in advance that, except for two 
extreme configurations, in general 0 H 0 M ; the relation between 0 H 
and 0 M has been discussed in detail in the preceding Sections). Note 
that the adoption of the formula (19) implies taking into account 
only two mechanisms of surface spin pinning: the isotropic pinning 
component a Q , the influence of which on the freedom of the spins is 
independent of their configuration with respect to the surface of the 
film, and the uniaxial factor a 2 (0 M ) representing the contribution of 
the uniaxial symmetry, with the surface normal as the symmetry axis, 
to the surface pinning. 

Already at this stage interesting conclusions regarding the prop- 
erties of the surface pinning can be drawn from the equation (19) 
despite its rather general formulation. Let us define two special 
angles: the critical angle 0 C M , for which natural pinning conditions 
prevail on the surface of the film, i.e., A(0) vf ) = 1, and the uniaxial 
pinning annihilation angle 0 M , for which the uniaxial pinning 
vanishes, i.e., 3 cos 2 9 M —1=0. The following equations apply to 
these special angles: 



a, 



(20) 
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A(6 u M ) = l-a 0 . (21) 

The latter equation provides a simple formula for the determination 
of the isotropic component a 0 of the surface pinning, only necessit- 
ating the value of the surface parameter in the external field config- 
uration corresponding to the uniaxial pinning annihilation angle 0 U M . 
With a 0 known, the configuration dependence of the uniaxial factor 
a 2 (0 M ) can be determined by the measurement of the surface para- 
meter A(0 M ) vs. 0 M (see the equation (19)). (We shall refer in this 
regard to the paper by Liu et al. 21 providing measurement data which 
will allow us to plot the experimental A(0 M ) dependence; see Section 
below.) On the other hand, theoretical considerations within the 
model used for describing the surface anisotropy in (Ga,Mn)As sam- 
ples will lead us to an equation, formulated in the next Section, in 
which a 2 (0 M ) is expressed by magnetic characteristics of the 
(Ga,Mn)As thin film; very interesting conclusions regarding the 
interrelation between the ranges of the exchange interaction on the 
surface and in the bulk of (Ga,Mn)As thin films will be drawn from 
the confrontation of the theory with the experiment. 



Confrontation of the Surface Pinning Model with 
SWR Study by Liu et al. 

Finally, the formula for the surface parameter takes the form: 

A(0 M ) = l-a 0 -[a 0 2 + a 2 {0 M )] (3cos 2 0 M -l), (27) 

where the coefficients a 2 and a\(9u) are as defined in (25) and (26). 
Note that in the surface inhomogeneity model the surface parameter 
(27) measures the degree of pinning of the surface spins and 
describes quantitatively the degree of the dynamic freedom with 
which they participate in the motion of the whole system of spins. 
The value A = 1 corresponds to a special case referred to as the 
natural freedom of the surface spins. Acquired by the surface spins 
as a result of breaking their interaction with those of their neighbors 
which are eliminated by the introduction of the surface, this freedom 
stems solely from the broken symmetry in the vicinity of the surface 
spins. Thus, absolute natural freedom of the surface spins only occurs 
when all the energy contributions in the equation (27) vanish simul- 
taneously, i.e.: 

a 0 = 0, (28) 



Model of the Uniaxial Surface Anisotropy in 
(Ga,Mn)As Thin Films 

We shall derive a phenomenological formula for the coefficient a 2 on 
the basis of our calculations presented above, in which the model of 
the uniaxial anisotropy is considered in both the microscopic and 
macroscopic approaches. From the equation (18) (see also Ref. [42]) 
it follows that the coefficient a 2 (0 M ) in the equation (19) can be 
expressed as: 



«2(0jh)= ^ 



4n AC 



-AC; 



surface 



J Dr- 



ill) 



where 47iM e «- = 47iM — H 2 ± , M is the saturation magnetization, 
H 2± the effective uniaxial anisotropy field, a the lattice constant (the 
average Mn-Mn distance), and D ex the exchange stiffness constant. 
The above equation indicates that both the intrinsic uniaxial aniso- 
tropy and the demagnetizing field contribute to the total uniaxial 
anisotropy in our model. 

As we will see later, extremely informative for the physical inter- 
pretation of the experiments performed by Liu et al. 11 is the express- 
ion of the latter contribution by the exchange length X, defined: 



hi 



Dex 



Dex 

4 n M sur f ac,! ' 



(23) 



we have introduced here a locally defined exchange length, different 
for the bulk and the surface. From the physical point of view it is 
reasonable to assume here that the lattice constant a in the equation 
(22) is identical with the exchange length ?. b that characterizes the 
interaction in the whole sample except for its surface. Under these 
assumptions (22) becomes: 



a 2 (0 M )=a° + a 1 2 (0 M ), 



\x\ 

2A,, 



a\(0 M ) = ■ 



Hi 



-surface 



-flff), 



(24) 



(25) 



(26) 



In the equations (24)-(26) we have indicated in advance what will 
follow from the confrontation of these formulas with the experi- 
mental data: that only the term a\(0M) is configuration-dependent. 



Hi 



-surface 



= H 



bulk 



Xb = x. 



(29) 
(30) 



However, as confirmed experimentally, the natural freedom of the 
surface spins is possible also in a particular situation in which the 
surface parameter value is one even though the above conditions are 
not all fulfilled. This particular situation may occur when there exists 
such a critical angle 0 L U that A(O^) = 1 because all the energy con- 
tributions in (27) annihilate each other. Further in this Section we 
shall analyze this situation in detail. 

On the basis of their SWR study of (Ga,Mn)As thin films Liu 
et al. 21 plotted the configuration dependence of the surface parameter 
A (0 H ) with the magnetic field rotating from the perpendicular {0 H = 
0) to parallel (f H = 90°) configuration (see Ref. [21]). As our formula 
(19) concerns the configuration dependence of the surface parameter 
versus 0 M , i.e. with rotating magnetization of the sample, the first 
thing necessary for proper interpretation of the measurements of Liu 
et al. is to find the dependence 0 M = 0 M (Q H ) in equilibrium condi- 
tions. The determination of the equilibrium conditions and the 
derivation of the sought relation 0 M = 0 M (0 H ) between the two 
configuration angles are presented above (see Figure 2). Figure 3 
shows the recalculated configuration dependence of the surface para- 
meter, with A plotted versus the new variable 0 M ; the plot corre- 
sponds to the measurement data of Liu et al. presented in Ref. 
[21]. The natural surface pinning is seen to occur for the critical angle 
6 M = 35° (which corresponds to the experimental angle 0 H = 19°). 
Also, the new plot reveals the occurrence of a local maximum in the 
A(0 M ) dependence around the angle 0 U M = 54.73°, for which the term 
(3 cos 2 0^ — l) equals zero (we shall take advantage of this finding 
below in further analysis of the experimental data of Liu et al. 21 .) 

Now we will demonstrate that the experimental curve shown in 
Figure 3 can be described by the function resulting from our SI 
model: 



A(6 M ) = l-a 0 -a 2 (0 M ) (3cos 2 0 M -l). 



(31) 



Knowing the maximal value of the surface parameter, 
A(6^A = 1.1068, we obtain immediately the value of the isotropic 
term in the series (31): 



a 0 = -0.1068. 



(32) 



On the other hand, the condition of occurrence of the local max- 
imum at 0 U U implies that the coefficient a 2 (0 M ) is zero at this point: 
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1= 0.85 

"> 20 30 40 50 60 70 80 90 

Magnetization angle, 0 M [deg] 

Figure 3 | Magnetization angle dependence of the surface pinning 
parameter A(0 M ) according to the experimental data obtained by Liu 
et al. 21 in their SWR study of a (Ga,Mn) As thin film; the plot corresponds 
to that shown in figure 9a in the cited paper, presenting the dependence 
on the magnetic field angle ff H . The applied transformation between the 
angles 6 H and 9 M is based on our determination of the equilibrium 
direction of the magnetization. 



« 2 (oy=o. 



(33) 



Both conditions allow to determine explicitly the function a 2 (Om) 
that reproduces the experimental plot shown in Figure 3 via the series 
(31). The determined function a 2 (0 M ) is presented in Figure 4. 

In the next step we shall refer to the formula (24) postulated in our 
model and representing the coefficient a 2 (0 M ) as the sum of a con- 
stant component a\ and a junction S 2 (0m). This implies that a 2 (Om) 
and 0^(0 m) have the same angular dependence, and their plots only 
differ by a shift a 2 along the ordinate axis. However, we do not know 
the value of a 2 \ This is a very sensitive point of our considerations, 
since in order to establish the value of a 2 we have to refer to the 
physical assumptions that are the very basis of our model of surface 
anisotropy. It seems reasonable to assume that of the three condi- 
tions (28, 29, 30) only (30) is fulfilled in the critical angle configura- 
tion; the other two energy contributions do not vanish, but 
compensate each other. This assumption means that by virtue of 
the equation (26) the coefficient a\ vanishes in the critical angle 
configuration: 

aj (<&)=<>, (34) 



which implies the equality: 



a\ = a 2 {0\ 



(35) 



Having established the value of the component a 2 we can already 
determine explicitly the function o^Ojh). The result is shown in 
Figure 5. 

From (26) it follows that: 



yJ\-2a\{0 M )' 



(36) 



and, on the basis of Figure 5, we can find the 0 M dependence of the XJ 
lb ratio. The obtained dependence is shown in Figure 6. Its analysis 
leads to very interesting physical conclusions. 

Note that in the plot in Figure 6 the surface exchange length l s is 
only slightly smaller than the bulk exchange length A;, for any angle 
0 M between the critical angle 6 C M and the parallel configuration angle 
0 M = 90°: 




0 10 20 30 40 50 60 70 80 90 
Magnetization angle, 6 M [deg] 

Figure 4 | Magnetization angle dependence of the surface pinning 
coefficient a 2 (ff M ) calculated from Eq. (31). 



0 C M <6 A 



<90 



(37) 



Thus, in this angle range a surface disturbance will not go beyond the 
first sub-surface plane formed by the spins directly under the surface. 
This means that the assumptions of the SI model are fulfilled very 
well in the angle range defined by (37)! In contrast, for angles 
0m < 0 C M X s i s greater than A b and grows steeply as the perpendicular 
configuration 0 M = 0 is approached. This means that in this angle 
range a surface disturbance, rather than being localized at the surface, 
penetrates into the bulk, affecting deeper sub-surface planes. Thus, 
the applicability of the SI model is very limited in this angle range, 
and the volume inhomogeneity model will be more adequate. This 
conclusion is fully confirmed by the experimental study by Liu 
etal. 21 . 

Further Physical Implications of the Model 

Now let us consider the component a 2 , which we have found to have 
a constant value, specified in (35), throughout the angle range 0 M e 
(0, 90°). The measurements performed by Liu et al. indicate that the 
material parameter values in the studied (Ga,Mn)As sample are D ex 
= 3.79 T • nm 2 and \nM eff = 4588 Oe, implying A b ~ 3 nm. On the 
other hand, for the critical angle Om = 0 C M from the formula (27) we 
get the equality: 
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Figure 5 | Magnetization angle dependence of the surface pinning 
coefficient ci\(Qm) calculated from Eq. (27). (See the text for detailed 
discussion.) 
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Figure 6 | Magnetization angle dependence of the IJl^ ratio resulting 
from our model of surface pinning in (Ga,Mn) As thin films (see Eq. (36)); 
2 S and 1^ denote the surface and bulk exchange length, respectively; 6 L M is 
the SWR critical angle. 



a 0 + a° 2 (3 cos 2 0 c M -l)=0, 



(38) 



which, after the substitution of a 0 = —0.1068 and 0 L M = 35°, yields the 
sought value: 

a°^0.108. (39) 

Now, getting back to (25), with the above-determined value of a\ we 
can estimate the difference between the effective uniaxial anisotropy 
field values on the surface and in the bulk: 



AH 2 



surface 
21 



H*f =913 Oe. 



(40) 



To our best knowledge, this is the first quantitative estimate of the 
surface uniaxial anisotropy field in (Ga,Mn)As thin films to be 
reported in the literature. 

As a measure of surface spin pinning experimentalists tend to use 
the surface anisotropy energy E s (0 M ), a phenomenological quantity 
thus related to the surface pinning parameter A used by us for 
describing the same feature: 

E s (0 M ) = M ^[A(0 M )-l]. (41) 

The above relation indicates that the character of the angular 
dependence of both quantities used for describing the surface pin- 
ning is identical, though in the equation (41) the reference level of the 
measure of the surface pinning is the zero value of the surface aniso- 
tropy energy, corresponding to our natural pinning A = 1. For E s 
(^m) > 0 the surface spins are unpinned, while for E s (0 M ) < 0 their 
freedom is constrained, which means that the surface spins are 
pinned. Plotted in Figure 7, E s (0 M ) has a maximum for Om = 0 M ; 
according to our estimate its maximal value is E S (0 M )~ 
0.07 erg/cm 2 . Note that this maximal surface anisotropy value is 
solely related to the free component a 0 , only responsible for the 
isotropic part of the surface spin pinning; the other surface aniso- 
tropy components only reduce this (maximal) value as the angle 
diverges from 0 " M in either direction. 

Final Remarks 

In our model the SWR critical angle is determined from the con- 
dition that the exchange length must be the same on the surface and 
in the bulk: 



Let us reflect on the physical mechanism as a result of which the 
magnetic properties of a (Ga,Mn) As sample fulfill the assumptions of 
the SI model for external field configuration angles 0h>0 c h , and 
those of the VI model for Oh < 0 C H . Let us get back to the formulas 
(23), which imply: 

(43) 



jy^surface 



M hulk 



This relation allows to determine, on the basis of the results presented 
in Figure 6, the configuration dependence of the surface to bulk 
magnetization ratio. The resulting dependence is shown in 
Figure 8. The plot indicates that if the sample is magnetized in direc- 
tions close to the surface normal (0 M ~ 0), the surface magnetization 
value is substantially reduced with respect to the bulk magnetization 
value. However, as the magnetization of the sample diverges from the 
surface normal to approach the orientation parallel to the surface 
(i.e., when 8 M ^> 90°), the surface magnetization (see Figure 8) gets 
closer to the bulk magnetization to equal it at the critical angle 0^, 
and surpass for larger angles, gaining an advantage stable throughout 
the angle range Om > 0 C M . The above interpretation is based on the 
implicit assumption that the magnetization of a (Ga,Mn)As thin film 
does not remain constant with changing out-of-plane direction 
(polar angle 0 M in the plane perpendicular to the surface of the film). 
Astonishingly, this assumption of ours proves fully consistent with 
the result of the very recently published Ref. [43]. In the cited paper 
the Authors show that the interpretation of their ferromagnetic res- 
onance measurements of bulk magnetization of (Ga,Mn)As thin 
films also requires the assumption that this magnetization is "a func- 
tion of the polar and azimuthal angles of the applied magnetic field", 
i.e., in our notation, M bu,k = M hu,k {0 H , <I> H ). In the context of such a 
result our above-presented considerations go beyond the conclusion 
of Ref. [43] by indicating that also the surface magnetization M'"^"" 
should be assumed to depend on the orientation of the external field; 
however, this dependence is different from that of the bulk magnet- 
ization. In the mechanism described above an important role should 
be assigned to the perpendicular surface anisotropy field acting along 
the surface normal. Note that the presence of this field also explains 
the existence of the surface mode (in the angle range 0 m > 0 C M ), since 
this very field gives the surface spins sufficient freedom for a spin- 
wave mode to be localized at the surface of the sample. 
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Figure 7 | Magnetization angle dependence of the surface anisotropy 
energy E s (fi M ) resulting from our theory (see Eq. (41)) for the (Ga,Mn)As 
thin film investigated by Liu et a/. 21 . 
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resulting from our model of surface pinning in (Ga,Mn)As thin films (see 
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Summary 

In this paper we show why in different conditions a resonance 
(Ga,Mn)As thin film sample may meet the assumptions of either 
the Surface Inhomogeneity (SI) or the Volume Inhomogeneity 
(VI) model. In our considerations we refer to the spin-wave res- 
onance (SWR) spectra measured by X. Liu et at 21 in (Ga,Mn)As thin 
films in different configurations of the static magnetic field H with 
respect to the surface. We demonstrated that the observed config- 
uration dependence of the SWR spectrum of the studied material can 
be described with the use of the surface pinning parameter expressed 
by the formula: 

A(0 M ) = l-a 0 -[fl° + ^(0 M )](3cos 2 0 M -l), 

where 0 M is the angle between the surface normal and the magnet- 
ization M of the sample. The values of the coefficients are estimated 
on the basis of the experimental data; the estimated value of the 
isotropic component of the surface pinning, a 0 = —0.1068, allows 
to determine the maximal surface anisotropy energy density, E s ~ 
0.07 erg/cm 2 . The intrinsic uniaxial anisotropy term a° 2 is of the 
order of 0.1, which implies that the uniaxial anisotropy field H 2 ± 
on the surface exceeds the bulk value by ca. 0. 1 T. We postulated that 
the coefficient a^Af) is related to the difference between the surface 
and bulk exchange lengths (A s and /if,, respectively). For a critical 
angle 0 C M , at which the SWR spectrum collapses to a single peak, X s 
= ).{,; the critical angle 0 C M separates two angle ranges in which the 
resonance properties are different: for Qm > 0 L M the SI model applies, 
since k s ~ /./,, and for 0m < 0 C M the VI model is adequate due to the 
domination of the surface exchange length (). s ^Xb). 
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